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We demonstrate that electroweak baryogenesis can occur in a supersymmetric model with an 
exact _R-symmetry. The minimal i?-symmetric supersymmetric model contains chiral superfields 
in the adjoint representation, giving Dirac gaugino masses, and an additional set of "_R-partner" 
Higgs superfields, giving .R-symmetric /i-terms. New superpotential couplings between the adjoints 
and the Higgs fields can simultaneously increase the strength of the electroweak phase transition 
and provide additional tree-level contributions to the lightest Higgs mass. Notably, no light stop 
is present in this framework, and in fact, we require both stops to be above a few TeV to provide 
sufficient radiative corrections to the lightest Higgs mass to bring it up to 125 GeV. Large CP- 
violating phases in the gaugino/higgsino sector allow us to match the baryon asymmetry of the 
Universe with no constraints from electric dipole moments due to f?-symmetry. We briefly discuss 
some of the more interesting phenomenology, particularly of the of the lightest CP-odd scalar. 



I. INTRODUCTION 

The origin of the matter asymmetry is a deep mystery 
that remains unsolved. Conditions that can lead to a dy- 
namical asymmetry between baryons and anti-baryons 
were articulated years ago by Sakharov pQ : baryon num- 
ber violation, C and CP violation, and a departure from 
thermal equilibrium. All three conditions are satisfied by 
the standard model as it passes through the electroweak 
phase transition. But, the CP violation is too small [2J, 
and the phase transition is not strongly first-order (e.g., 

Weak scale supersymmetry has long been known to 
potentially enhance the strength of the electroweak phase 
transition and provide new sources of CP violation [8HTU] . 
In the minimal supersymmetric standard model (MSSM), 
this occurs in the presence of a light stop and a light 
Higgs boson. Given the recent LHC results interpreted 
as the existence of a Higgs boson at rrih = 125 GeV [TTJ 
fl2] , this region is essentially ruled out [T3l [14] . Methods 
to strengthen the first-order phase transition beyond the 
MSSM have been widely discussed [T5H3I]- 

In this paper we consider a relatively recent framework 
for supersymmetry that incorporates an i?-symmctry, 
proposed in Ref. |35j . /("-symmetric supersymmetry fea- 
tures Dirac gaugino masses, that have been considered 
long ago |36H38| and have inspired more recent model 
building [351 EHHHS] and phenomenology (3TJ I51HT5] . We 
show that an i?-symmetric supersymmetric model can si- 
multaneously obtain: a strong enough first order phase 
transition; sufficient CP violation with no difficulties with 
electric dipole moment (EDM) bounds; and a Higgs mass 
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m/j ~ 125 GeV consistent with the LHC observations. 
Much of these results rely on exploiting the additional 
superpotential couplings among the Higgs fields, their 
i?-symmetric partners, and the chiral adjoint fields. Ku- 
mar and Ponton also studied electroweak baryogenesis in 
a model with an approximate i?-symmetry [31] . Their 
approach was to reshuffle the i?-charges of the fields such 
that ^gH u H^ operator is allowed, where the fermion sin- 
glet in $b is the i?-partner to the bino. In our approach, 
we retain the original i?-charges defined by the minimal 
i?-symmetric supersymmetric standard model [35 , utiliz- 
ing new superpotential couplings among the electroweak 
adjoints, the Higgs superfields, and the i?-partner Higgs 
superfields. 

Supersoft supersymmetry breaking [39] shares several 
ingredients of the i?-symmetric model. One positive fea- 
ture is the relative weakness of the all hadronic jets plus 
missing energy search bounds from LHC due to the heavy 
Dirac gluino mass |75| . On the flip-side, however, the 
usual D-term that determines the tree-level contribution 
to the lightest Higgs mass is absent, and no A-terms are 
generated. Thus, even with some nontrivial modifica- 
tions to the model [59] , it seems rather difficult to rec- 
oncile the recent LHC observations of = 125 GeV 
[TT1 fT2"] with the predictions of the supersoft model. In 
contrast, one of central points of our paper is to show that 
there are tree-level (and loop-level) contributions to the 
Higgs mass from the same superpotential couplings that 
allow the electroweak phase transition to be strength- 
ened. These additional contributions imply i?-symmetry 
need not be broken to generate a large enough lightest 
Higgs mass. However, we will still need a substantial one- 
loop contribution from stops with mass ~ 3 TeV to ob- 
tain TTih ~ 125 GeV, and so some sacrifice in fine-tuning 
is inevitable. 



II. THE MINIMAL R-SYMMETRIC 
SUPERSYMMETRIC STANDARD MODEL 

First we review the field content and new couplings 
present in the minimal i?-symmetric supersymmetric 
standard model (MRSSM). In the MRSSM, the gaugi- 
nos acquire Dirac masses through the Lagrangian terms 
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where W™ is the field strength superfield for one of the 
SM gauge groups (labelled by a, a is a spinor index) and 
$ a is a "i?-partner" chiral superfield transforming under 
the adjoint representation of the appropriate gauge group 
with i?-charge i?[$ a ] = 0. Supersymmetry breaking is 
communicated through i?-symmetry preserving spurions 
that include W' a which parameterizes a -D-type spurion, 
W' a =D8 a . Expanded into components, the above op- 
erator becomes 
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(X a ^ a + h.c. + V2D a {A a + A a *)) = 
Id (A ^ a + h.c. + 2 V2 D a Re(^)) , (2) 



that contains the mass term between the gaugino (A a ) 
and its "i?-partner" (tp a ) as well as a coupling of the real 
part of the scalar field within $ a to the D-term of the 
corresponding gauge group. 

The second term in Eq. ([2| has two important con- 
sequences: First, the equation of motion for Re(A a ) sets 
D a = for all three SM gauge groups. The Higgs quartic 
coupling in the MSSM is contained in the SU{2) and U(l) 
D-terms, so eliminating these terms will clearly have an 
impact on the Higgs potential. Second, while the real 
parts of A a acquire a mass 0{M D ) from Eq. lm{A a ) 
remains massless at this level. 

In order to enforce i?-symmetry on the superpotential, 
the Higgs sector of the MRSSM must be enlarged. The 
/i-term of the MSSM is replaced by the i?-symmetric /i- 
terms 
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where R u ,d are new, i?-charge R[R U: d] = 2 fields that 
transform as (1, 2)^/2 under the standard model gauge 
groups. This choice of i?-partners ensures that elec- 
troweak symmetry breaking by the Higgs fields H u d 
does not spontaneously break i?-symmetry. The MRSSM 
also defines the i?-charges of the matter fields to be 
R[Qi, U[, Df, Li, Ef] = 1, allowing the usual Yukawa cou- 
plings in the superpotential. 

Given the extra matter content, there are new super- 
potential operators |35) one can write in the i?-symmetric 
theory, 

W D \ u B <S> B H u R u + \ d B <$> B R d Hd 

+ \^ a w H u T a R u + \ d w <S> a w R d T a H d . (4) 



Unlike the ^-terms, which are required to achieve exper- 
imentally viable chargino masses, there is no direct phe- 
nomenology that dictates that the A^ couplings in Eq. Q 
must be nonzero (being superpotential couplings, they 
will not be generated radiatively if set to zero initially). 
However, these couplings play a vital important role 
in driving the phase transition to be first order. The im- 
portance of the Xi couplings can be seen already from the 
scalar potential; the operators in Eqs. (3j4| lead to new 
trilinear and quartic operators involving Higgs fields and 
the scalars in <I> B , & w , R u ,R d . 

V D ^(\: B A* B )\Hy + ^ d (K B A B )\H° d \ 2 +c.c. : (5) 

Trilnear scalar interactions involving the Higgs multi- 
plets, especially those with large couplings, are well 
known to impact the strength of the electroweak phase 
transition [TU1 HI 1221 HE [S3 [MSI] • 

Turning to the supersymmetry breaking parameters of 
the theory, scalar soft masses can arise from an additional 
source of _F-term supersymmetry breaking. So long as 
the supersymmetry breaking spurions X have i?-charge 
R[X] — 2, the i?-symmetry is preserved and no Majorana 
gaugino masses are generatedjj The soft masses from the 
Kahler terms are 
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In addition, holomorphic soft masses for each $ a are of 
the form 
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We assume the coefficients for the holomorphic soft 
masses are real. The full set of soft masses for the scalar 
components of $b and <E»vk ar e given in the Appendix 
in Eq. ( A5 1 . Soft-breaking, trilinear scalar couplings be- 



tween the Higgs and squarks or sleptons are forbidden 
by i?-symmetry. For viable phenomenology, we allow the 
relative size of the supersymmetry breaking contributions 
to be within roughly one order of magnitude in mass. 

Throughout this paper we will take the Dirac gaug- 
ino masses to be large. This limit simplifies our calcu- 
lations and is motivated by phenomenology. Specifically, 
to avoid conflict with precision electroweak observables 
Ref. [35] found the SU(2) W gaugino masses should be 
larger than 1 TeV. Such heavy electroweak gauginos de- 
couple from the rest of the theory and play little role in 
the electroweak phase transition. The higgsino masses in 
the MRSSM, on the other hand, come from \i Ul /id, which 
we take to be closer to the electroweak scale. 

Furthermore, heavy Dirac gauginos, when combined 
with the MRSSM Higgs superpotential structure and 



1 i?-symmetry is not essential here. Majorana gaugino masses can 
be avoided as long as X is not a singlet |50| . 
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lack of A-terms, lead to significantly relaxed flavor con- 
straints. As shown in Ref. |35) . low-energy, precision ob- 
servables become insensitive to new sources of flavor or 
CP in the supersymmetric sector. Electric dipole mo- 
ments induced by one-loop contributions involving the 
gauginos and higgsinos are completely absent. This al- 
lows 0(1) phases in the MRSSM that will be important 
when we consider CP violation and its role in baryogen- 
esis in Sec. IVII 

Having reviewed the MRSSM, its typical spectra and 
constraints, we now investigate the strength of the elec- 
troweak phase transition. 

III. QUALITATIVE FEATURES OF THE PHASE 
TRANSITION IN THE MRSSM 

There are several features of the MRSSM which lead 
one to suspect that the phase transition could be different 
from more familiar (supersymmetric) scenarios. First, as 
a result of the superpotential interaction in Eq. Q, there 
are extra scalar states coupled to the Higgs boson. Extra 
"higgsphilic" scalars are known to (potentially) increase 
the strength of the EW phase transition, with the proto- 
typical example being the stop squark. However, unlike 
the stops of the MSSM, these MRSSM states are not col- 
ored, and they have limited interactions with other SM 
fields. As a result, these extra scalars can be quite light 
and can interact strongly with the Higgs without causing 
any phenomenological problems. 

The second feature is that the tree-level Higgs poten- 
tial vanishes to leading order in 1/Md where Mjj is the 
Dirac gaugino mass for the bino and wino. This can be 
understood as follows: in the limit that all other super- 
potential or soft-breaking interactions involving Re(A a ) 
are absent or negligible, the Z?-term disappears. In the 
MSSM, the D-term is the sole source of Higgs quartic 
interactions at tree-leveQ D — means the potential 
is purely quadratic and tree-level symmetry breakdown 
does not occur. In the presence of other interactions in- 
volving Re(A a ), the Z?-terms are not exactly zero and a 
(non-trivial) tree-level Higgs potential is generated. The 
dimension of the Higgs operators that are generated de- 
pends on how Re(A a ) interacts, but all operators will be 
accompanied by coefficients with at least one power of the 
large Dirac gaugino mass, 1/Md. As a simple example 
demonstrating this mechanism, consider the potential 

V(a,h) = ml h 2 + (M D a - gh 2 ) 2 + \\iah 2 + m% a 2 . 

(8) 

Though simpler than the full MRSSM potential, this toy 
potential has all the important features; a Higgs boson 



h and a scalar a that couples to a "D-term" g h 2 . We 
can study this potential in three interesting limits: If 
A/i = mg = 0, the field a can be integrated out exactly 
and the residual potential is purely quadratic in h. If 
A/i 7^ 0, ms = 0, an effective Higgs quartic is generated 

Notice that the quartic receives a positive contribution (if 
A/i > 0) to order 0(1/Md), and a negative contribution 
to order 0(1/M D ). We will see the same result in the 
MRSSM, suggesting a modest hierarchy with [i <C M D 
maximizes the quartic coupling. If ms ^ 0, A/i = (while 
assuming ms <C Md), we get a different effective Higgs 
quartic 

A.W(f|) + c(4). (10) 

Finally, if both A/i ^ 0, ms ^ 0, the effective quartic is 
the sum of the last two equations. We can remove the 
quadratic term by demanding a minima at h = v; the 
resulting potential is then entirely proportional to A e //, 
and therefore oc l/M^^n > 1. Explicitly, 

V = ml h 2 + X eff h 4 — > X eff {h 2 - v 2 ) 2 , 

where we have added an unimportant overall constant. 

Because the tree-level potential is suppressed, the 
nature of the zero temperature electroweak symmetry- 
breaking (EWSB) minima of the full 1-loop potential is 
somewhat more complicated. The EWSB minima can be 
nearly degenerate, or even higher than the zero temper- 
ature, trivial vacuaj^] When the trivial vacuum is only 
slightly higher than the EWSB vacua (at T = 0), the 
critical temperature will be low, making large 4> C /T C eas- 
ier to achieve. A cartoon depicting the nearly degenerate 
minima scenario with more typical scenarios is shown be- 
low in Fig. [T] This effect is mitigated somewhat by the 
presence of an effective Higgs potential below the scale 
of the Dirac gaugino masses (needed to obtain a phe- 
nomcnologically viable Higgs mass). 



IV. EFFECTIVE POTENTIAL IN THE MRSSM 

In this section we will describe the various scalar field 
contributions from the MRSSM that enter the effective 
potential at both zero temperature and finite tempera- 
ture. 



2 Strictly speaking, this assumes only renormalizable superpoten- 
tial terms are included. Higher dimensional operators will change 
this statement, as studied in Ref. 82-85 . 



3 The EWSB conditions used to fix the soft masses m 2 Hu ,m 
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only ensure the electroweak vacua is a local minima, and not 
necessarily the global minima. 
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FIG. 1. Cartoon showing the T — Higgs potential in a 
nearly-degenerate minima scenario, as in the MRSSM (black 
solid) and a more conventional Higgs potential (solid green). 
The zero of these T — has been chosen such that V(cf> c ) = 
0, where, roughly, (f> c is the Higgs vacuum expectation value 
(vev) at the critical tempterature. When the temperature is 
raised slightly, the potential shifts to the dashed lines. For 
the degenerate case, the small shift in temperature brings us 
to T c , while for the typical potential the temperature must 
be raised higher for a phase transition. A high temperature, 
T 3> all masses, potential is also shown (blue). 



A. Zero temperature Potential 

We will make the following simplifications to avoid an 
overabundance of parameters: (Dirac) gaugino masses 
Mi = M2 = M, " fi -terms" fi u — /id = A*, additional 
Yukawa couplings Af = Af = As, Ajf = A^ = Xw 
and equal soft-masses for R u ,d 1 mR u — m Rd = m R- To 
satisfy precision electroweak constraints, we take M to 
be much larger than all other scales. The neutral scalar 
components of $s and §w ar e expanded as 
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The first ingredient in the calculation is the tree- 
level potential Vt ree , which consists of the superpotential 
piece, the soft masses and the D-term potential. As we 
discussed earlier, the interplay between the D-term and 
the other contributions results in a viable Higgs potential 
with an EWSB vacuum. 

Faced with the hierarchy among superparticle masses 
(M 3> jU, etc.), we proceed by integrating out all particles 
with mass M. These include the gauginos (both charged 
and neutral) and several scalars. The scalars with mass 
M include two CP-even neutral states and one charged 
scalar. The origin of their mass can be traced back to 
Eq. pj): they are, up to small mixing effects, the real 
scalars within As and A l w . Removing the heaviest fields, 



the residual potential now has 0{l/M) and 0(l/M 2 ) 
suppressed interactions. 

The full potential, and the resulting tree-level mass 
matrices, is shown in Appendix [Xj The mass matrices in 
the low-energy effective theory are kept field-dependent, 
meaning we retain all h u = (f> u , hd = 4>d dependence. 

Focusing on the CP-even, neutral scalar sector, we next 
calculate the 1-loop Coleman- Weinberg (CW) potential. 
Working in the MS scheme, 
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where 8V is a counterterm that is added to Vcw to en- 
sure that the 1-loop potential has an extrema at <p u = v u , 
4>d = v d- We have made the distinction between <p u d 
and v U} d since, at finite temperature, the fields will de- 
viate their T = vacuum values (it is also convenient 
to parameterize the Higgs field dependence in terms of 
4> 2 = (0?, + <^d) aim X = arctan((/) It /0 d )). The sum in 



Eq. ( 12 ) runs over the relevant (light) particles, with m 



counting the degrees of freedom, m {4> u ^(j)d) representing 
the field-dependent mass, and Ci is a constant equal to 
3/2 for fermions and scalars and 5/6 for gauge bosons. 
All calculations have been performed in Landau gauge 
(for recent discussion on gauge-dependent artifacts, see 
Ref. [86]). 

We will come to the exact states included in the sum 
shortly, however note that the field content is slightly 
different than in the MSSM. The sum in the MRSSM 
contains no gauginos, but includes all i?-partner fields 
(both the scalar and fermionic components). To sim- 
plify the calculation, we will neglect sleptons, first and 
second generation squarks, and the sbottoms since their 
couplings to the Higgs are small. The final parameter 
in the CW potential is the renormalization scale A. In 
order to minimize the effects from higher order terms, we 
take A equal to the mass of the heaviest dynamical field, 
A = max(ffli(^)) [57] , 

Before moving to finite temperature, the total (tree + 
1-loop) T = potential must satisfy several consistency 
checks. First, the EWSB minima must be the lowest 
minima in order for the vacuum to be stable rather than 
meta-stable. 



V T=0 ( 



, 4>d = v d ) < V t =q{4>u =4>d = 0) (13) 



While this condition is always applied, it has little impact 
on the parameter space of models with an unsuppressed 
tree-level potential. The second condition is that the 
EWSB is a minimum and not a saddle point. The coun- 
terterms added to Vcw only require the vacuum values 
0ii = v Ul (j)d = Vd extremize the potential. To assure a 
minimum, we must also enforce 
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This condition is automatically satisfied so long as the 
Higgs boson masses are positive. 

The effects of the various fields on the Higgs poten- 
tial clearly depends on their mass and spin. Under the 
assumption that M 3> all other mass scales, the mass 
eigenvalues fall into several categories: 

Category 1: The first category contains particles of 
mass ~ M ; very heavy fields that we have already inte- 
grated out. 

Category 2: The second category contains lighter 
fields that have 0-independent mass. These fields shift 
the Higgs potential only by an overall constant and are 
therefore unimportant to our calculation of the phase 
transition. The higgsinos fall into this category, as do 
a full multiplet of Higgs scalars (charged, neutral CP- 
even, neutral CP-odd). The Higgs multiplet in this cate- 
gory corresponds roughly the H°, A and H ± of the usual 
MSSM. Because of our assumption that the R u .d scalars 
have a common mass, one of the neutral scalars also has 
a ^-independent mass. 

Category 3: The third category contains fields with 
mass of the form m + /(</>), where m is a weak scale 
parameter (/i, or one of the soft masses other than M) 
while the Higgs field-dependence is an additive func- 
tion f{4>). The remaining charged Higgs scalars, the 
imaginary parts of the Ab,A 1 w scalars, the charged R- 
Higgs scalars, and one of the neutral i?-Higgs scalars 
all have masses of this type. Being light and with <j)- 
dependent masses, these states are especially relevant for 
us, so their masses are explicitly displayed below using 
<$>u = 4> sinx, <j>d = <f> cosx- 
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m| + A 2 4> 2 sin 2 x 
= H 2 + m 2 R + A 2 4> 2 cos 2 x 
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where the soft masses 
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m 2 p0l m 



m pt for the scalars 



in Eq. (Ill are defined in the Appendix in Eq. (A5). 

Category 4: The fourth and final category contains 
fields whose mass comes entirely from electroweak sym- 
metry breaking. This includes the weak gauge bosons 
and the top quark. If certain combinations of the Ag, A l w 
soft masses (i.e. A a) happen to be small, one or more 



of Hf,A%, A% will also receive their mass entirely from 
electroweak symmetry breaking. 



B. Lightest Higgs Mass 

The lightest Higgs boson mass in the MRSSM deserves 
special attention. In general, it receives three main con- 
tributions to its mass to one-loop: 
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(16) 



Unlike the MSSM, there are no tree-level contribu- 
tions from the usual Z?-term [33]. This would-be dis- 
aster is averted in the MRSSM due to new tree-level 
contributions from the A-terms as well as soft-mass 
contributions to the adjoint scalars. The general ex- 
pression can be straightforwardly evaluated numerically 
from the effective potential, which we do in our nu- 
merical results below. The leading contributions, to 
C(l/M 3 ), can be obtained analytically in the limits 

4/ 



\n/M\, v/\M\, m 2 /M 2 < 1, and tan/3 > 1: 



m 



fa. tree 



V2M 

v 2 
32M 2 



Mz/j.(X\v cos9w — 2As sin 8w) 



4A|(3Mf sin 2 W - 4/z 2 ) 

- A 2 v (4 M 2 + 3M|cos 2 6 W ) 

- &X B X w M 2 z sin 2 V 



sin 2 9w(m 2 + m 2 pQ ) + cos 2 d w {m 2 t + m 2 t ) 



pi i 



4M 2 



O 



AP J 



(17) 



The tree-level contribution is maximized when Xwf 1 > 
simultaneously with As/i < 0, with a phase convention 
where the Dirac gaugino masses are real and positive. 

Taking A = X w = -As, m s0 = m p0 = m st = m pt = 
ms, and evaluating the contributions for characteristic 
masses that we will see later in our numerical evaluation: 

2 n71f 2i A 4 1T eV 

fc.tree Z 2 Q0 GeV M 



~ §.\%M 2 7 X 2 



200 GeV 



2 / 1 TeV 



M 



0.02Af 2 A 2 ' 1TcV V 



IT) 



l -M 2 ^- 
2 Z A/ 2 



(18) 



This approximate expression slightly underestimates the 
tree-level contribution. Nevertheless, we see that we can 
achieve m\ tree nearly equal to M§ with A ~ 2, intrigu- 
ing (though accidentially) similar to the largest tree-level 
value found in the MSSM. 

The one-loop contributions to the lightest Higgs mass 
from the stops <5m 2 - are well-known [55] and won't be 
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repeated here. We note, however, that there is no scalar 
trilinear coupling A t tj,t* R h due to i?-symmetry, and thus 
At = 0. This means we will need stops with masses above 
a few TeV to obtain a large enough one-loop radiative 
correction to the Higgs mass, and as a consequence, we 
can integrate the stops out in the calculation of the elec- 
troweak phase transition. 

Finally, there are additional one-loop contributions 
to the Higgs mass from the terms proportional to A. 
These are straightforwardly derived from the Coleman- 
Weinberg potential. 



C. Finite Temperature Contributions 

The effective potential at T > can be separated 
into a temperature-independent contribution as well as a 
temperature-dependent contriubtion. The temperature- 
independent contribution is the tree-level plus Coleman- 
Weinberg potential calculated in the previous section. 
The temperature-dependent contribution includes 



J±(y) 



rp4 

n i^2 J i( m i(<l>u,<l>d)/T 2 ), 



(19) 
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where J+(J_) is the thermal function for fermions 
(bosons), respectively. The thermal potential must be 
amended due to some well-known subtleties of pertur- 
bation theory in finite temperature, however before ad- 
dressing these it is important to break down the effects 
of moving to T ^ 0. 

In the limits T 3> m and T <^ m the thermal functions 
have a simple form 



V T = 



no 



\2tvT 2 J 



3/2 



2 / T T < m 



(20) 



In light of these limits, finite temperature effects from 
particles with mass M T are completely negligible. 
Similarly, fields with purely ^-dependent mass have the 
largest impact on the shape of the potential. For = 0, 
these fields are light (up to thermal masses, which we will 
come to shortly), so the thermal contribution is negative 
definite and oc T 4 , meanwhile, out at larger field values 
(of 4>) , all these fields are heavy so the thermal corrections 
are small. Thus, fields whose mass comes entirely from 
<j) push the trivial (<f> = 0) vacuum downwards sharply as 
the temperature increases while leaving the large-0 part 
of the potential unaffected. Because the thermal correc- 
tions at <t> — depend so strongly on T, the more degrees 
of freedom in the m oc category, the lower we need to 
raise the temperature before the trivial vev and EWSB 
vev equilibrate, leading to larger </> c /T c . For fields with 
mass of the form m = juq + f((f>) the thermal contribu- 
tion depends in detail on the relation between m and T, 
so we must evaluate these contributions numerically. 



1. Thermal masses 

Thermal masses are systematically calculated by sum- 
ming over the "daisy" diagrams [5 : , 89, 90 where the con- 
tribution to their mass is typically of order g 2 T 2 . Phys- 
ically, they represent the screening of scalar fields in a 
thermal bath. Their effect is to reduce the ^-dependence 
in scalar field masses, hence weakening the phase transi- 
tion. 

Thermal masses are most important for fields whose 
mass is determined entirely from electroweak symmetry 
breaking (m oc cf>), since these fields become massless as 
(\> — > 0. For fields whose mass has a (/(-independent piece, 
the thermal contribution has little effect. Therefore, we 
include thermal masses only for the charged Higgs, JJg 
and the relevant CP-odd charge- neutral Higgs fields A® 4 . 
Thermal masses for the longitudinal W/Z degrees of free- 
dom (0(g 2 )) and all other scalars (in categories 1-3) are 
neglected for simplicity. 

We evaluate the thermal mass correction in the interac- 
tion basis. We use the expressions in Ref. [HI], ignoring 
terms proportional to the electroweak gauge couplings 
and only retaining terms of order A 2 T 2 , with A ~ O(l). 
In the large gaugino mass limit, the only fields we need 
to consider are the Im(As) and the In^A^); the light 
pseudoscalars are combinations of Iib.(Ab) and lm(Ay V ), 
and the charged Higgs fields are made up of Im(^4^,). For 
this subset of fields, the thermal mass enters as additional 
terms in the interaction basis mass matrix Hij, where i, j 
runs over (Im(A ), Im(A 3 )) ((lm(A+), Im(vl~*)) + c.c.) 
for the CP-odd Higgs (charged Higgs). Only self en- 
ergy diagrams giving rise to a quadratic divergence at 
T = contribute to the thermal mass (i.e. requires 1/p 2 
in the integrand of the loop integral by power count- 
ing). Therefore, only self energy diagrams originating 
from the 4-point scalar interactions and scalar-fermion- 
fermion interactions without interior chirality flips con- 
tribute. By this argument, supersoft interactions do not 
contribute to the thermal masses as they contain only 
3-point scalar interactions. Similarly, Z?-term contri- 
butions are proportional to the electroweak gauge cou- 
plings and they are negligible in comparison to A. The 
only interaction terms that can contribute come from 
Eq. Q. However, the interactions in Eq. ^ do not 
generate off-diagonal thermal masses - the correction to 
the Im(Ao) — Im(^43) entry is proportional to = 0, 

and there no sufficiently strong/divergent interactions to 
create a Im(A + )lm(A~'*) entry. The diagonal elements 
Haoa°,^-a 3 a 3 ,^a+a~ ^a+*a~' i on the other hand, are 
non-zero: 



(21) 
(22) 

A^ T 2 (23) 



Ha°a° 


= -\ 2 T 2 
2 B 


n A 3 A 3 


1 \2 rp2 

- r wI 


Ha+a- 


= Ha+*a- 



G 



Iii the mass basis, 



Eq. (16 1, removing all (^-independent contributions to 



\2 

-2 2 i A W t>2 



(24) 
(25) 



To incorporate the thermal masses, we follow the pro- 
cedure in Ref [511551150] and include a "ring" contribution 



V r i n g = — V rii (mf (<f> Ui <j> A ) - rhi (</> u , cf> d )), 

iZTT z — ' 

i 

where i runs over {A®, A° A , H 3 }. 



seal. 



(26) 



V. PHASE TRANSITION: 
RESULTS 



NUMERICAL 



Over the last few sections, we have discussed the com- 
ponents of the complete, finite temperature plus 1-loop 
potential 



V = V tr 



V cw + V T + V, 



ring • 



(27) 



We are now ready to assemble the pieces and be- 
gin our hunt for regions where the phase transition is 
strong. Starting with the tree-level Higgs potential, we 
integrate out all heavy fields; for the mass ~ M scalars 
and gauginos, this is done at tree level, while the stop 
squarks must be integrated out at one-loop level. This 
potential, restricted to light fields, is then augmented 
by the Coleman- Weinberg potential. We only include 
states in the Coleman- Weinberg that are light and whose 
couplings to the Higgs are unsuppressed - specifically, 
SM gauge bosons, the top quark, and the six states in 
Eq. (16). The CW contributions from </>— independent 



states merely shift the effective potential by an additive 
constant and are uninteresting. For a given parameter 
set, the scale A in the CW potential is set to maxfm^)) 
to minimize the effect of higher order corrections. The 
finite temperature potential - for the same set of states 
as we used in the CW piece - is then added, along with 
ring contributions for any bosonic fields whose mass is 
directly proportional to </> (i.e. gauge bosons). 
As in Sec. 



IV 



we will make several assumptions in 
order to reduce the parameter space for our numerical 
study: 



Mi = Mi = M 
Hu = = A* 



\u \d \ \u 

A R — A R = A R , A w 



\;; A/; = A/;. A(( — A'/ ( = \\\ 



(28) 



m Ru 



L s0 



4 



: m Rd = m R 

2 _ _2 _ 
l p0 



The first five assumptions are identical to those in 
Sec. |IV| The final assumption, the equality of the F- 



term masses for the A« scalars sets A,, 



A'. 



and m 2 A o (in terms of the mass categories laid out 

4 

earlier, these two states move from category three to cat- 
egory four). Under these assumptions, the strength of 
the phase transition is a function of eight parameters: 
tan/3, rn,A,M,fx, Xb, Xw,fns anc ^ m R- However, since 
mji always appears with /1 in the combination /i 2 + m R , 
we will set = and use /j as a proxy for the two, 
reducing the problem to seven parameters. 

Requiring the gaugino mass M > 1 TeV, we scan over 
the other parameters, evaluating 4> C ,T C and the light- 
est Higgs mass at every point. Starting with the full 
Higgs potential raised to a sufficiently high tempera- 
ture (^250 GeV), we lower the temperature in succes- 
sively smaller steps until we find a minima in the po- 
tential that is degenerate with the electroweak symmet- 
ric minima. In principle we ought to search for min- 
ima in a two-dimensional space: (<f> u , 4>d) OI \ equivalcntly 
((ft, tan%). However, in performing this two-dimensional 
search, we find that the ratio of the Higgs fields at T c 
remains close to the vacuum value (the ratio of vevs) 
provided tan /3 is > 3. Therefore, to more quickly scan 
over several parameters, we focus on tan /? > 3, set 
tanx = tan/3(T c ) = tan/3 = tan/3(T = 0), effectively 
reducing the Higgs potential to a function of cj) alone. 

With M and tan/3 essentially fixed by our approxi- 
mations, our parameter set is reduced to m^,/x, As, Aw 
and m 2 s . Of these, and m 2 s have little impact on 
the strength of the phase transition; both appear only in 
the tree-level Higgs potential (given our assumptions in 
Eq. ( 29 )), either as (^-independent terms or suppressed by 
two powers of the heavy scale M. Therefore, we show our 
results by fixing itla — 300 GeV, ms — and scanning 
in various directions over the most sensitive parameters, 
fi, Xb, Xw- 

Our first scan, shown in Fig. [2j shows the dependence 
on the Xi couplings for fixed /i. This particular scan was 
performed using M = 1 TeV, tan /3 = 4, though increas- 
ing to larger values for either parameter has a negligi- 
ble impact. Within the (Ab,A^) space, the 4> C /T C con- 
tours trace out a "bullseye" shape. This is expected - 
the largest effect on the strength of the phase transition 
should come when the interactions between the Higgs 
fields and new scalars are strongest. The difference in 
normalization between Xb and Xw, along with factors of 
tan 9, set how fast the contours change in Xb versus Xw- 

The second panel of Fig. [2] shows the lightest Higgs 
mass in the same parameter space. Clearly, the region of 
Xb < 0, Xw > is favored since it allows rrih > 125 GeV. 
The shape of the Higgs contours is driven by the Aj de- 
pendence of the tree- level piece, which we explored in 



Sec. IV A Loop level contributions to mi, though siz 



in 



able, do not prefer a given sign for the A^ since they are 
always proportional to X B ,X W . We emphasize that the 
relative sign of Xb,Xw is, of course, convention depen- 
dent. Specifically, the signs of superpotential couplings 
depend on the ordering of the fields in Eq. ([3j [4]) . 

In Fig.[3j we zoom in on the most interesting quadrant, 
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FIG. 2. On the left, we show contours of C /T C in the (X B ,Xw) plane with fj, = 200 GeV and M = 1 TeV. The other 
parameters chosen were tan/3 = 4, = 300 GeV, m 3 = GeV, but the values of <j> c /T c are not particularly sensitive to 
these choices. On the right, we show contours of the lightest Higgs mass (in GeV) assuming m^ L = m,^ = 3 TeV. Large stop 
masses were necessary to obtain contours that approach ran = 125 GeV. 



Xb < 0,\w > 0, overlaying the lightest Higgs mass con- 
tours on top of the 4>c/T c contours. To demonstrate the 
effect of the Dirac gaugino mass, we repeat this zoomed- 
in scan for a second Dirac mass, M = 2 TeV. For larger 
M, the strength of the phase transition is hardly changed, 
while the mass of the lightest Higgs is slightly reduced 
since the tree-level contribution to mf L scales as 1/M. In 
both panels of Fig. [3] we can see that there are regions 
where the Higgs mass is close to 125 GeV and the phase 
transition is strong, (f> c /T c > 1. 

To study how the strength of the phase transition de- 
pends on [a, we consider a different direction in param- 
eter space, where A^ = — A\y while we scan over Ayy 
and /i. The resulting <j) c /T c and contours (overlaid) 
are shown in Fig. [4] For the parameter ranges we have 
plotted, the /i dependence of the phase transition is mi- 
nor. As in Fig. [3j we see that there are regions where 
a strongly first order electroweak phase transition can 
be achieved simultaneously with a 125 GeV Higgs boson. 
The viable parameter regions all require A ~ 1.5 — 2. 
Had we extended the scan in Fig. [4] to larger values, we 
would find regions where the condition in Eq. ( 14 ) is vio- 



lated because the lightest Higgs mass (squaredjls driven 
negative. 

Summarizing our numerical studies, we have shown 
that the electroweak phase transition can be strong over 
a wide range of viable MRSSM parameters. However, a 
strongly first order phase transition is only part of the 
baryogenesis mechanism. A first order phase transition 



ensures baryon-number-violating sphaleron process are 
out of equilibrium. This prevents sphalerons from eras- 
ing any generated baryon asymmetry, but we still need 
to generate an asymmetry in the first place. One well- 
established way to generate an asymmetry is through CP- 
violating collisions between particles in the plasma and 
the walls of vacuum bubbles. We explore this mechanism 
in the context of the MRSSM in the next section. 



VI. CP-VIOLATION 

Tunneling processes can be understood semiclassically 
by spacetime-dependent field configurations that connect 
the real and false vacua. For the EW phase transition, 
the field that interpolates is the Higgs field, 4>(x). Ex- 
panding about the tunneling configuration, fields (top 
quarks, charginos, etc.) that interact with the Higgs ap- 
pear to have a spacetime-dependent masses; if the Higgs 
is complex, these masses will also be complex. Com- 
plex masses violate CP, so collisions of particles with this 
spacetime-dependent Higgs can be shown to lead to an 
asymmetry between particles and antiparticles on either 
side of the wall [52Tl95j . A CP asymmetry generated 
at the outer edge of the bubble wall is then propagated 
to the interior of the bubble (where sphaleron effects 
are unsuppressed) via diffusion and transport mechan- 
ics [HJ 15514100] . Provided a strong enough source of CP 
violation, the observed baryon asymmetry can easily be 
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FIG. 3. Same as upper-left part of Fig. [2] except we have "zoomed in" to the contour regions of interest. The left-side plot has 
M — 1 TeV, as in Fig. [2] while the right-side plot has M = 2 TeV; all of the other parameters are the same as given in Fig. [2] 
Here we have overlaid Higgs mass (in GeV) contours to demonstrate the region of overlap for mj, ~ 125 GeV simultaneous 
with <j> c /T € > 1. 




FIG. 4. Contours of <f> c /T c (left panel) 4> c (in GeV, middle panel), T c (in GeV, right panel) in the (A, /i) plane. In these figures 
we have taken A = X w = -As and ^ = 200 GeV, M = 1 TeV, tan/3 = 4, m A = 300 GeV, and m s = GeV. 



generated. The difficulty lies in introducing a CP violat- 
ing source that is strong enough, once manipulated into 
the form of a complex Higgs coupling/complex mass, yet 
secluded enough to avoid conflict with existing flavor and 
CP observables. Along with the magnitude of the CP- 
violating phase, the thickness and speed of the bubble 
are important parameters in determining the size of the 
asymmetry. 

The most studied MSSM baryogenesis model uses 
the relative phase between the /x-term and the gaug- 
ino mass (A/2) as the CP- violation source [10] . How- 
ever, this phase is stringently bounded by EDM measure- 
ments [101] . and successful MSSM baryogenesis relies on 



resonance production when /i ~ M^. In the MRSSM, we 
have the opposite problem. There is no significant EDM 
constraint on the CP phase due to the absence of Ma- 
jorana mass and left-right squark mass mixing |35 [ I54 j . 
but there is also no resonance production enhancement 
as M2 ^ [i is required by EW precision measurements. 
A more careful study is therefore necessary to show that 
the MRSSM can generate enough asymmetry. 

There are several complex parameters in MRSSM: two 
higgsino mass terms /i u and /i^; three Dirac gaugino 
masses M^; three holomorphic scalar masses of the ad- 
joints mf ; the term; and four A couplings in Eq. Q, 
which give 13 complex phases. The phases of seven su- 



perfields H u ^, Ru,d, &b w g can remove srx of the phases 
while keeping an i?-symmetry one. The seven surviving 
phases can be parametrized ag^j 

m t M*,i = B,W,g 

Mu Mj (XiY , Md Mj (A J d ) * , j = B, W. (29) 

Under our simplifying assumptions of equal A", X d cou- 
plings, equal electroweak gaugino masses and equal 
[ii terms, the number of phases is reduced to three: 
Arg(m s Mt),Arg([j,M X B ),Arg(nM X w ). The strong 
CP constraint sets a bound Arg (m^M-) <C 10~ 7 , but 
there are no significant bounds on the /Lt-term related 
phases. 

The phases Arg([x M Xb), Arg([i M Xw) can source CP 
violation in the MRSSM through the higgsino sector. At 
leading order in l/M, the higgsino masses are simply /i. 
However, after integrating out the mass M gauginos and 
scalars, the higgsino masses shift to 



g ( /) 2 (z)e 2tq * ( - z '> 



16M 



- sin 2 (3 (2X B tan 6 + X 



w ) 



+ cos 2 /3 (2X B t&nO + (1 + 4V2)X X 



(30) 



g<\> 2 (z)e 2lq ^ 



16M 



cos 2 /3(2A B tan6> + X w ) 



sin 2 (3 (2X3 tan 6 + (1 + 4v / 2)A w ) 



where we have added z dependence to quantities that 
will vary across the bubble walJ^J and subsumed the 
relative phases [iMXw, [iMXb into a single (spacetime- 
dependent) phase for <f){z). The phase in the higgsino 
masses is suppressed by l/M, however, as emphasized 
above, we can consider much larger phases (and phase 
changes) Aq ~ tt in the MRSSM as there is no EDM 
constraint. 

To get a better estimate of the asymmetry that can 
be generated from chargino interactions, we rely on the 
similarity between the higgsino masses in the MRSSM 
and the higgsino masses in the nMSSM. Specifically, as 
shown in Ref. [21], in the limit Mi 3> /x, the dominant 
source of CP violation in the nMSSM comes from the 
changing phase of /1. Within this limit and assuming 
canonical profiles [B] for the phase and magnitude of /_t 
across the bubble, Ref. [H] found the generated baryon 
to photon ratio to be: 



7/10 ^ c(T c ) 



Aq 1 



TT r . 



.3/2 



A/i 



exp 



(31) 



4 See Appendix [X] for the relation between the soft masses 
m a t,m a Q etc., and holomorphic/non-holomorphic soft masses. 

5 For simplicity we will assume a planar bubble propagating in the 
z direction. 



where 7710 = 10 10 77, c(T c ) ~ 1.6T c /GeV, r ~ 0.78, fx is 
the amplitude of the z-independent part of the higgsino 
mass, and A/j (Aq) is the change in amplitude (phase) 
of the coordinate-dependent part of /i across the bubble. 
The remaining parameters in Eq. (31) are the thickness 
of the wall l w and the critical temperature T c . 



Applying Eq. (31 1 to the MRSSM, we make the fol- 



lowing identifications 



Ho = M 



(32) 



A/j 
Aq 



iv 2 cos 2(3 



AM 



(2X B tan9 + (l + 2V2)X w ), 



(33) 



where Afj, is derived by summing over the coordinate 
dependent parts of both chargino masses. The A/j, also 
includes an extra factor of two to account for the fact 
that the /j D 4> 2 {z) in the MRSSM, while /1 D <j)(z) in the 
nMSSM [21] , thus the rate of change across the bubble 
is twice is large. 

As an illustration, consider the point Xw = — = 2, 
M = 1 TeV, fi = 200 GeV, m A = 300 GeV, tan /3 = 4, 
and m s = GeV. The critcal temperature can be read off 
from Fig.[4]to be T c ~ 135 GeV. Using these parameters, 
including only the chargino contributions to the baryon 
asymmetry, with l w = 10/T C [21] and Aq = tt, we obtain 
i]i Q ~ 4.0. This is larger than needed to match the baryon 
asymmetry of the Universe, but it trivial to adjust Aq 
or other parameters to bring this in line with 7710 ~ 1. 
We have not considered possible contributions from Dirac 
neutralinos, nor from squarks and sleptons, which would 
be interesting to study in future work. 



VII. COLLIDER LIMITS 

As expected from phase transition lore and shown 
explicitly in Sec. |IV C[ the phase transition becomes 
increasingly first order when light scalars are present. 
Therefore, in order to judge how strong the phase tran- 
sition can actually be in the MRSSM, we need to know 
just how light the relevant scalars can be. We begin with 
a recap of the light particles in this scenario (full mass 
matrices and eigenstates can be found in Appendix |A|) : 

• the lightest neutral Higgs boson (other Higgs 
scalars have mass ~ 777^4) 

• two CP-odd scalars, linear combinations of the 
pseudoscalars in the W3 and B i?-partners: we la- 
bel these A%,A% 



• one of the charged Higgs scalars, H. 

• higgsinos and i?-higgsinos 



Five of these six states reside at the weak scale and 
are light only by comparison to the gauginos. For the 
(R-)higgsinos, charged Higgs, and one combination of 
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CP-odd scalars {A®), this weak-scale mass comes at tree 
level, ~ [i for the fermions and ~ A v for the scalars. 
(Note that for the range of A we are interested in, this is 
high enough to avoid any direct bounds.) As we showed 
earlier, the Higgs boson can be made sufficiently heavy 
through the combination of a small 0(Xfi/M) tree-level 
contribution and large radiative corrections. The mass of 
the remaining light state A3, another CP-odd scalar, is 
insensitive to the Higgs vev and is instead set by a differ- 
ence in soft masses. Because this state is independent of 
the Higgs vev, it does not play a direct role in setting the 
strength of the phase transition, so its phenomenology 
may seem unrelated to the issues in this paper. However, 
the same combination of soft masses that enters into m A o 

is also present in the mass of the and A° - two states 
that play a large role in strengthening the phase transi- 
tion. To increase the strength of the phase transition, 
H ± ,A\ should be as light as possible, which means we 
want to take the soft mass contribution to m H ± , m A o to 

be small. For this configuration of soft masses, A® will be 
much lighter than the weak scale (and possibly massless 
at tree level) , so the viability of this parameter set given 
collider bounds is far from guaranteed. 

The A® is a combination of gaugino i?-partner fields, 
so it does not interact directly with SM fermions. The 
wino i?-partner does have gauge interactions, while the 
bino iZ-partner does not. However, we are only interested 
in the neutral wino (T3 = 0) i?-partner component - 
this component does not interact with at all and 
only has a quartic interaction with the charged SU(2) 
gauge bosons. Furthermore, the mass eigenstate (A3) 
coupling is suppressed by a mixing angle (squared) ~ 1/2 
as only the wino i?-partner component has this quartic 
interaction. 

The remaining interactions are: a small trilinear 
W^W~ ,1 (A°) term, suppressed by ~ v/M, and Yukawa 
interactions with higgsinos and i?-higgsinos. These in- 
teractions play no role in A3 production, but they do 
permit A3 to decay through loops. The only final state 
that can proceed through a higgsino loop is A3 1 — > 77. 
The gauge boson loop can lead to photons as well as 
light fermions, but the fermion component is subdomi- 
nant due to the suppression by Yukawa couplings (and 
the overall coupling ~ v/M). For simplicity, we therefore 
assume A3 — > 77 is the only available decay mode. 

Given its interactions with standard model particles, 
it is quite challenging to produce the A®. Because the 
A3 does not interact with the Z , LEP places no con- 
straint. The only production mechanism (at leading 
order in v/M) is vector boson fusion (VBF), e + e~ — > 
vv + W + W~ — > vv + (A]) 2 , leading to the spectacular 
final state of missing energy plus four photons. However, 
even assuming a branching ratio to photons of 100% for 
A", the rate is much too small to have been seen at LEP. 

At hadron colliders, there are two ways to produce A3 1 ; 
through VBF: pp,pp — > jj + (A3) 2 — > jj + 47, and from 
s-channel W production followed by the emission of a 
pair of A3 1 : pp,pp -> W ± + (A3 1 ) 2 -> W ± + Aj. Provided 



some of the photons have high-py, such signals would be 
clearly visible above background. The issue is whether 
the rate is high enough to generate more than a handful 
of events. 

At the Tevatron, the rates (at leading order) for jj + 
(A§) 2 before VBF cuts are - 0.5 fb for m A o, 3 = 20GeV, 
falling to - 0.2 fb for m A o 3 = 50 GeV. At the LHC 
(8TeV), we find a(pp -> jj{A°) 2 ) LO ~ 26, 19 fb for 
m^o ,3 = 20, 50 GeV respectively. The significant increase 
in the rate at LHC is because qq can initiate VBF pro- 
duction. After standard object identification and fiducial 
volume cuts (not even the usual VBF cuts) and account- 
ing for realistic identification efficiencies (even if we look 
for fewer objects, like 77+ jj), the Tevatron rates are too 
low to provide any bound. We arrive at the same conclu- 
sion for the W ± + (A3) 2 process; the rate at the Tevatron, 
while slightly higher than the VBF process (~ 0.7 fb for 
tti^o 3 = 20 GeV) is still too low to provide a meaning- 
ful bound given the Tevatron dataset and realistic object 
efficiencies. 

At the LHC, the VBF rate is high enough that a more 
thorough investigation is necessary. Di-jet plus multi- 
photon events would certainly fall under the scrutiny 
of the LHC Higgs searches. To test bounds on A§ 
coming from Higgs diphoton limits, we generate signal 
events using the machinery of MadGraph4 |102j . and 
Delphes |103| . then pass events through a mock CMS 
Higgs analysis. Though CMS and ATLAS have ded- 
icated "VBF" searches looking for 2 jets and 2 pho- 
tons [104| I105j , the analysis looking for a final state most 
similar to the jj(A§) 2 final state, we find the VBF cuts 
imposed are too restrictive and hence the signal efficiency 
is extremely low < 10~ 4 . The more inclusive diphoton 
Higgs searches have looser cuts, but we also find them to 
be not particularly sensitive to our signal, e ~ 10~ 3 . The 
lack of sensitivity is due to a few reasons: While there are 
more photons in our signal, the photons themselves have 
lower energies and sit in a more crowded environment 
as opposed to the diphoton signal (~ 125 GeV Higgs) 
the cuts were designed for. Hence the leading photon pt 
cut, photon-jet isolation, and photon-photon isolation re- 
quirements remove much more signal compared to a SM 
Higgs. A heavier A3 would pass the cuts more efficiently, 
but has a smaller production cross section. 

While we find no firm bounds on the A3 from Higgs (or 
other) searches, the LHC rate is certainly large enough 
that a dedicated multi-photon plus jets (or plus W^) 
may well be worthwhile. 

VIII. DISCUSSION 

We have seen that electroweak baryogenesis can be 
achieved in the minimal supersymmetric i?-symmctric 
model with: 

• an electroweak phase transition strength 4> C /T C > 1 

• Higgs mass ~ 125 GeV 
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• induced baryon asymmetry rjio > 1. 

The central ingredients are the new superpotential cou- 
plings, Eq. Q, where we required Xw — —X B ~ 2 to 
achieve a strong enough first order phase transition si- 
multaneous with to/j ~ 125 GeV. 

That we needed modestly large As providing substan- 
tial trilinear interactions between the Higgs boson and 
the additional scalars in R u , R d , $s and &w is perhaps 
not particularly surprising, given the degree of freedom 
counting given in Ref. [T7]. Larger A couplings are po- 
tentially problematic if the theory is run to higher scales, 
though this is beyond the scope of this paper. However, 
there are a few comments we can make on this point. In- 
terestingly, the interactions between the chiral adjoints 
$ a and the Higgs/i? superfields are also present in mod- 
els with Af = 2 supersymmetry involving TV = 2 vector 
supermultiplcts interacting with Af = 2 hypermultiplets 
(e.g. 1106] ). There the coupling strength is deter- 
mined by Af = 2 supersymmetry to be \[2g for the ap- 
propriate gauge group (times for fj(l)y). This is 
somewhat smaller than, but not that far from the su- 
perpotential coupling strengths of interest in our case. 



If we had taken the Af = 2 limit for the superpoten- 
tial couplings, they would evolve identically to the gauge 
couplings up to the explicit Af =2 breaking terms. This 
suggests that the renormalization group evolution is not 
necessarily as drastic as, say, the superpotential coupling 
for SH u H d in the NMSSM. It would be interesting to 
investigate this further, and to determine the role of the 
relative signs of these couplings on the evolution. 
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Appendix A: Scalar potential 

The contributions to the scalar potential comes from the superpotential, the D-term, supersoft terms and the scalar 
soft masses. The superpotential is 



W RSSM = {Hi R-) (/i„ + X uB A - ^A 3 ) - (H° u R° u ) (/z u + X uB A + ^A, 



X 



(R+ H~) ( Md + X dB A Q - '-^A 3 ) - (R d fl£) ( fM + X dB A 



A- 



A~ ( \ A + / 

-j= [X uW (Hi Rl) + X dW (R d H° d )j — [X uW (H° R~) + X dW (R d H~) 



c.c. 



(Al) 



From the superpotential, we get the usual potential V D — |Fi| 2 , where i runs over all of the above superfields. 
The D-term contribution to the Lagrangian is 



D-term 



I if. 



+ ' 2 + |tfT 



|i?;| 2 -K| 2 "l^°| 2 -|^| 2 



lie 



l^l 2 ) 5 



2 

H^ H° u + R°*R- + H-*H° d + R+*R° d - 2iA* 2 A 3 + c.c 



iH Q u *Hi + iR-*R° u + iH^*H° d + iR°/R+ - 2iA%A l + c.c. 



(A2) 
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The supersoft terms are 

supersoft = -M 2 (A + A* ) 2 V2g'M 1 (A + A* )h\H+\ 2 + |tf°| 2 \H d \ 2 \R~\ 2 - \R° U \ 2 + |i?+| 2 + \R° d \ 2 ) 



-M 2 ^ + A*) 2 - V2gM 2 (A 3 + A* 3 )±(\H+\ 2 - \H°J 2 + \H° d \ 2 - \H d \ 2 \R^\ 2 + \R° U \ 2 
-^2gM 2 {Ai + Al) l -{H+*H° U + H°/H d + R° U *R^ + R? ' R° d + c.c.) 
-V2gM 2 (A 2 + A*)^(-iH+*H° - iH° d *H d - iR° u *R- - iR+*R° d + c.c.) 



We further decompose the neutral fields A and A3, CP even and CP-odd pieces 



\Rl\ 2 -\R°, 



(A3) 



. _ sp+ipo _ s 3 + ipa 

A °-^T' Ai -^r- (A4) 

The CP even pieces will mix with the CP-even Higgs scalars h u ,h d , while the CP-odd pieces only mix among 
themselves. 

The soft mass terms are the usual MSSM soft masses, plus equivalent R u ,R d soft masses. The only tricky soft 
masses are for the Aj 

2 2 
V soft D m 2 (A Q A*) + m f(A 2 + Ay 2 ) + m 2 t (AA*) + ™f(A 2 + A*' 2 ), (A5) 



where i is an SU(2) index. In the language of Sec. VI m 2 ,rn 2 t are holomorphic soft masses (and can potentially 



carry a phase), while m 2 ,m 2 t are non- holomorphic, and therefore purely real, soft masses. 

Appendix B: Light field potential 

Starting with the full potential given in the previous appendix, we first integrate out the scalars with mass M, 
keeping terms of O(l/Mo) and 0(1/M^), The Higgs scalar soft masses m 2 ^, rn 2 Hd can be removed by enforcing 
electroweak symmetry breaking. Specifically, under the set of assumed relations between various parameters laid out 



in Eq. ([29]), we find. 

,2 



2 2 22^ uMzv sin 8(cos6\w — 2sin0A.s) ./ 1 \ 

l H + M = m A COS P 7= ^ & "T-77T 

Hu ah 2V2M \M 2 ) 

2 2 22^ uMzv cos 2 (3(cos9 \w ~ 2sui0Ab) / 1 \ 
lHd +H =m A cos p ^- + °{W 2 )- 



(Bl) 



Though we have only shown the modifications to 0(1/M), we retain terms to 0(1/M 2 ) when calculating the tree-level 
potential. Focusing on neutral, CP even Higgs fields and plugging in the expressions for m 2 H , m 2 Hd , we find the tree 
level potential to be 

T/ m \ ,2 .2/0 \ 1 l*Mz 4> 2 {{4> 2 -« 2 )cos2x- « 2 cos2/3) . 

Vtree sin (p — x) H 8V2Mv (2As sin 6 — Xw cos 6) 



(2M 2 Z cos 2/3 cos 26> (cos 2/3 + 2cos2x)(4A| - \ 2 W + k\ B \w tan20) 



1024M 2 

-(4A| + \ 2 W ){M 2 Z (1 +4 cos 2(3 cos 2 X + cos 4/3) -32 ^ 2 ) 

-64M| cos 2/3 cos 2x (sin 2 6>(to 2 + m 2 ) + cos 2 6»(m 2 t + m 2 t ))) 

(2A/| cos 2 2 X (sin 2 9{m 2 pQ + mg ) + cos 2 8(m 2 pf + m 2 t )) - ^ 2 v 2 (4A| + A^)) 
^ 64M 2 u 2 

6 M| cos 2 2x(2A B sin 6 - \ w cos 6>) 2 
^ 256 M 2 u 2 ' 



(B2) 
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where (f) u , 4>d are the up- and down- type Higgs fields and (f> 2 = <$>\ + (j) 2 ,, tan% = 4) u /<Pd 

In addition to the tree-level potential, our calculation also requires the field dependent masses for all of the light 
states (mass <C M) in the theory. By field-dependent we mean that the Higgs fields cj) u and 4>d arc not set to the 
zero-temperature vacuum values v u ,Vd- The mass matrices and eigenstates are given in the subsequent subsections. 
In all expressions we neglect any 0(1/M) or smaller pieces. This truncation is justified because the effects of these 
states on the Higgs potential are already suppressed by loop factors. 



CP-odd, Charge Neutral Higgs Scalars 



The field-dependent mass matrix for the four R = 0, neutral, CP-odd fields (a u , ad,po,P3) is block diagonal and is 
given below: 

/ m A cos 2 /3 m 2 A cos j3 sin /3 \ 





m 2 A cos (3 sin (3 m 2 A sin 1 (3 





l s0 








^ + <fi) m 2 t -m 2 pt + x -f- {4> 2 U + <p 2 d )) 



(B3) 

The pseudoscalar Higgs mass matrix has a zero eigenvalue corresponding to the Goldstone boson eaten by the Z. 
The three massive eigenvalues are 



m 2 A o 2 


= m A 






m A°,3 


= A A + 




4- 1 X 2 
~r 4 A W 




4 


m 2 A o A 


= A A + 


X 2 


-1- 1 X 2 
+ 1 A W 




4 



{<t>l + 4>l) - \[*!% + A' A (A| - \\ 2 W )(4>1 + 4> 2 ) + (A| + \\ 2 W )^ 2 U + 4> 2 d )\ 
{4>l + 4>l) + J A' 2 + A' A (\ 2 B \\ 2 W )(4>1 + 4>l) + (A| + \x 2 w W u + <!> 2 d ) 2 , 



A A = m 2 s0 + m 



A A = m 



l p0 



(B4) 



up to corrections of order 1/M. 



Charged Higgs Scalars (R = 0) 

For the charged Higgs scalars, we started with four states that all mix with each other: ff+, Hj* , A + , A*~ + c.c. 
One combination of A + , A*~ is heavy and gets integrated out. The remaining three- by-three mass matrix is 

m 2 A cos f3 2 m A cos j3 sin f3 
m 2 A cos/3sin/3 m 2 A sin/3 2 

m 2 st - 

This mass matrix has a zero eigenvalue corresponding to a Goldstone degree of freedom. The remaining eigenvalues 
are 

m H±, 2 = m A 

m 2 H±3 =m 2 st -m%+ > f( ( f>l + ( l> 2 d ) (B6) 




Neutral R — 2 scalars 



The mass matrix for the complex, neutral scalar i?-partners of the Higgs fields is: 

/ M 2 + m 2 ? + A|+ j A ^ e/> 2 u \\ B Xw^u4>d 2 \ 
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The eigenvalues become particularly simple if we take m Ru = m R d = m-n 

X 2 -4- IX 2 

2,2 2 i 2 i ~ A W ri2 , ±2\ 

m R + fi , m R + fi + ^ + 

in which case only one of the states has a </>- dependent mass. 



(B7) 



Charged R = 2 scalars 

The charged i?-scalars do not mix with any other states, so they have a simple mass term 



( f + < + & <&) \R-f + (m 2 + m R + >f 4>l) |i?+| 2 



(B8) 



CP even neutral scalars 



Though we start with four CP-even, charge neutral scalars (h u , hd, sq, S3), two have mass ~ M and are integrated 
out of the low-energy theory. The remaining two states to form the light Higgs boson h and its heavier cousin H . 
The heavy Higgs boson has mass ~ tjia and plays no role in determining the strength of the phase transition. The 
light Higgs boson is massless at lowest order at O(M ), but is lifted to nonzero mass (at tree level) by 0(1/ M, 1/M 2 ) 
terms in Vt ree (Eq. (B2 1 ). The full expression is long and not very insightful. The mass matrix simplifies significantly 



at large tan/3, leading to the expression given in Eq. (jTTJ) . 

Appendix C: Ino mass matrices 



The neutralino mass matrix starts as 



M 2 








Mi 






w/ 2 

A i 4>d 




2\f2 


V2 



_gjh L g<t>d \ 

, 2 2 

a 4>u a >Pd 

2 2 



-Mi- 








(CI) 



Under our parameter assumptions, once the mass M scalars and gauginos are integrated out, the matrix collapses to 
a diagonal two-by-two matrix with entries ±/z + 0(1/ M). 

The chargino mass matrix combining W + ,RJ with A~,HJ is originally 

Sri 

After our usual parameter assumptions and upon integrating out the heavy states, we are left with a single entry 

M z 



< 2 ^ 

X d Vd 

2 



m 2 ± 

Xi 



H - 



\M 



- (2A B sin 0(tf d - 2 ) + cos 6((1 + 4^2)0 2 - 2 ))) 



(C2) 



where we have retained the 0(1/ M) piece within m 2 ± because it will be relevant for the CP-violation section. The 



Xi 



extra 2\/20 2 in the above expression comes from integrating out the gauginos (at tree level), in addition to the heavy 
scalars. 

Similarly, the matrix combining A + ,H^ with W~,R~ begins as: 



M 2 
V2 



2 

/■'.« 



(C3) 
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becoming 



m 2 ± 
x 2 



M z 
8Mv 



2Xb sin( 



4>l) + cos 8{4> 2 d -(1+W2)^)) 



(C4) 



once the heavy states are removed. 
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